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SUMMARY 


, . * fv,,;, fifteen elastic constants associated 

Formulas are derived for the shear w of plates with closely 

with bending^ stretching, v s | configurations and proportions, 

spaced integral ribbing in a variety of configuration ^ ^ 

In the derivation the plates ara ^°^ der of’plywood. The constants, which 
orthotropic plates somewhat on resisting deformations other 

include the effectiveness of f ? lo^iSdinal^irections, are defined 
than bending and str ®^ h ng n „ , ^ and theoretical and 

in terms of four coefficients a, P, * , ^ s /^ eff lcients are dis- 

experimental methods for e e i as tic constants are predicted by 

cussed. Four of the more Import ^elastic const correlation is 

these formulas and are compared with test resuxus. 
obtained. 


introduction 


. , Qt 1 f f ened construction, evidenced by 

Growing interest in integr y large forging press program 

such papers as references ^n J or production, emphasizes the 

i£-fS S C o h ^“on P rS characteristics of integrally 

stiffened plates. 

, . +(a for the prediction of structural characteristics 

A primary requisite fo p ^stic constants. In the present 

of plates is a knowledge of t fifteen elastic constants 

paper, therefore, formulas ^sti^ shearing of pistes 

associated with the bending ^ more directions. The 

with closely spaced integral formulas are illustrated in figure 1 and 
ribbing patterns covered by the fomulas a 3ect ion is arbi- 

section rib Sth circular fillets at its base. 
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The elastic -constant formulas derived involve four coefficients a, 
S, a' , (3' for each rib which define the effectiveness of the rib in 

resisting deformations other than simple bending or stretching in its 
longitudinal direction. For most purposes a reasonably accurate evalua- 
tion of these coefficients is required. Experimental and theoretical 
methods of evaluating them are discussed. 

As a check on the correctness of the elastic-constant formulas, the 
predictions of the formulas for four of the more important elastic con- 
stants are compared with experimental data. 

The principal symbols used are defined in appendix A. 


DEFINITION OF ELASTIC CONSTANTS 


If the rib spacings are small compared to the plate width and length, 
it is plausible, for purposes of studying over-all or average behavior, 
to assume that the actual plate may be replaced by an equivalent uniform 
orthotropic plate. Figure 2 shows an infinitesimal element of the equiv- 
alent plate subjected to bending moments of intensity Mx and My, 
twisting moments of intensity M^, stretching forces of intensity N x 
and Ny acting in planes I and II, respectively, and shearing forces of 

intensity N^y in plane III. The locations of planes I, II, and III are 
arbitrary. 

The behavior of the element can be described by a set of force - 
distortion relationships in which elastic constants appear. Such rela- 
tionships are obtainable from reference 5* If deflections due to depth- 
wise shear are assumed to be negligible as is customary in ordinary plate 
theory, the following equations (eqs. (l') to (6') of ref. 5) are obtained 
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where -— and — - are the curvatures, ^ is the twist, e x 

dx dy 

and €y are the extensional strains in planes I and II, respectively, 
and /xy is the shear strain in plane III. 


According to these equations, fifteen constants are needed to estab- 
lish the force-distortion relationships - namely, two bending stiffnesses 
D x and D y , a twisting stiffness Dxy, two stretching moduli E x and 
Ey, a shearing modulus Gxy, two Poisson's ratios p x and My associated 
with bending, two Poisson's ratios p' x and p' y associated with 
stretching, four coupling terms Cxx, Cxy, Cyx, and Cyy associated 
with bending and stretching, and one coupling term T associated with 
twisting and shear. Not all these constants are independent, however, 
for,, as a consequence of the reciprocity theorem for elastic structures, 

Py = %Px/ D x M-'y = % P'x/ E x* 

The form in which the force -distortion relationships have just been 
given is not the most convenient form for some applications, for example, 
for buckling calculations. For such purposes a more suitable form is 
obtained when the first three equations are solved simultaneously for 
M x , My, and Mxy and these expressions are then used to eliminate M x , 

My, and M^ in the last three equations. The six new force-distortion 

equations thus obtained are 
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M xy " 2D k + C k N xy 
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where Hy = I>2 M-x/^l an d ^2 = ^2^l/^l* 

Of the fifteen elastic constants appearing in equations (7) to (12), 
two, n x and n , were also in the original set of force-distortion 

equations. The remaining constants (D^, D2, Dk> ®i, E£, Gk> nq, 

1^2, C 1;L , Cj_ 2 ’ C 21 , ^ 22 ) 811(1 Ck) are' new. The algebraic relation- 
ships between the new and the original elastic constants are given in 
appendix B. 


METHOD OF ANALYSIS 


The analysis is made for a plate with the general pattern of ribbing 
shown in figure 3(a), which includes, as special cases, the patterns of 
figure 1. A typical repeating element of the plate is indicated by the 
short-dashed rectangle in figure 3( a ) and is shown three-dimensionally 
in figure 3(b) • 

The analysis is based on the assumption that each of the four rib 
segments shown in figure 3(b) may be replaced by three orthotropic sheets 
of material parallel to the skin, each one covering the entire area b x by 

and each fastened to the skin by means of many hypothetical, perfectly 
rigid, infinitesimally small bars imbedded perpendicularly through the 
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bending in its longitudinal direction, another (labeled (2) ) represent 
onS^he effectiveness of the rib in resisting stretching and bending in 
its transverse direction, and the third (labeled © ) represents only 
the effectiveness of the rib in resisting shearing and twisting relative 
to its longitudinal and transverse directions. (The transverse direc- 
tion, as used herein, is the direction in which t w is measured see 
fig. 5.) In order for the three substitute sheets to accomplish their 
purpose, they are assigned the following properties: 

(a) Sheet © has a volume equal to that of the rib segment it 
replaces, with its center of gravity at the same level as that of the 
rib Its stretching or compressing modulus of elasticity in the direc- 
tion of the rib is E and its modulus transverse to the rib is zero. 

Its stiffness per unit width for bending in the direction of the rib 

is equal to the bending stiffness of the rib about lts 

py the rib spacing (i.e., b x for a y-wise rib, by for an x-wise rib, 
and b s for a skew rib, fig. 3(a)), while its bending stiffness in the 
direction transverse to the rib is zero. The shearing and twisting 
stiffnesses and Poisson's ratios of the sheet are assumed to be ze . 

(b) Sheet © has a volume equal to some fraction p of the volume 
of the rib segment, with its center of gravity at some distance aH 
above the middle surface of the skin. The modulus of elasticity for 
stretching or compressing in the direction transverse to the rib is 
whereas that in the longitudinal direction of the rib is zero. The bending, 
shearing, and twisting stiffnesses, and Poisson s ratios for sheet (2) 

are 8.11 assumed "to he zero. 

(c) Sheet © has a volume equal to some fraction p' of the volume 
of the rib segment, with its center of gravity at some distance a H 
above the middle surface of the skin. Its modulus of alas *J^ y f °£ h 
shearing relative to the longitudinal and transverse directions of the 
rib is G, whereas its twisting stiffness relative to these two dire c- ^ 
tions is zero, as are the stretching and bending stiffnesses and Poisson s 

ratios . 

It is assumed that the hypothetical sheets offer no interference to 
one another. 

On the basis of the foregoing assumptions, the integrally stiffened 
plate has been converted to a more homogeneous plate somewhat on the 
order of plywood. The assumption of rigid bars connecting the substitute 
sheets and the skin is equivalent to the assumption that material lines 
normal to the surface of the plate before deformation remain straight 
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fjnyl Tig rj^f nrTnaf.l on . If it is f UT Lher assumed that these lines remain 
perpendicular to the surface of the plate and that the stresses are in. 
the elastic range, any of the methods used for ordinary isotropic plate 
analysis may he readily extended to the present idealized structure. 

For the present purpose an energy method is adopted to determine 
the six forces and moments necessary to maintain the prescribed uniform 

deformations e x> e y> ^ > ami 7xy’ equations 

Sx^ by x 

obtained for these forces and moments in terms of the distortions are 
put in the form of equations (l) to ( 6 ) to yield formulas for the 
original elastic constants or in the form of equations ( 7 ) to ( 12 ) to 
yield formulas for the new elastic constants. 

The details of the analysis and the derivation of the elastic 
constants are presented in appendix C. The formulas obtained for these 
constants are presented in the following section and the evaluation of 
a, p, a', and p' is discussed in two succeeding sections. 


FORMULAS FOR ELASTIC CONSTANTS 


In this section the formulas are presented for the calculation of 
the fifteen elastic constants appearing in equations 1 to 6 and the 
thirteen new constants appearing in equations 7 to 12. The formulas 
are presented for the most general type of plate considered, which is 
illustrated in figure 3 . For plates with one or more sets of ribs 
omitted, the formulas also apply when the terms representing the areas 
and moments of inertia of the omitted ribs are set equal to zero. 

The formulas for the constants in the original force-distortion 
equations (l) to ( 6 ) are as follows: 
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1 

2 ( k xy - k m) 

EH 2 

I xy _ 


(27) 


The formulas for the constants in the new equations (eqs. (7) to (12)) 
are as follows: 
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where 

b x , by, b s the spacing of the x-wise, y-wise, and skew ribs, respec- 
tively, in. 

9 the angle of skew of the ribbing, deg 

H the over-all height of skin plus ribs, in. 

tg the thickness of the skin, in. 

p Poisson's ratio for material 

Equations (43) to (54) contain the quantities A Wx , A Wy , and A Wg , 
k Wx , k Wy , k Ws , and I Wx , I Wy , and I Wg which define the areas, loca- 
tions of centroids, and moments of inertia of the rihs. For rectangular 
ribs with circular fillets, as shown in figure 4 , these quantities are 
given by the equations 
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( 56 ) 


( 57 ) 


(Eq. (57) contains a factor 2 to account for the fact that there are two 
ribs in the skewed direction - one at an angle +6 to the x-direction 
the other at an angle -9 to the x-direction.) 


= 


W 1 



/ WxA h/ 


tW x t S + l^S m 
2 H 


t S V 


H 



16 



NACA RM L53E13a 


The values of kj, kjp and kjjj depend upon the locations of 
the centroids of the forces N x , N y , and Nxy, respectively, imposed 
upon the plate element. (See fig. 2.) For the important case in which 

N x acts in such a plane that it produces no curvature 


and 


acts in such a plane that it produces no curvature 
must equal zero (see eqs. 1 and 2) and, therefore. 


d 2 w 


dx 
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Similarly, for the case in which Njjy acts in such a plane that it 

d 2 * 

produces no twist , T must equal zero and, therefore, 

dx dy 


'III - k xy 


(66) 


If N x and N y do act in such planes that they produce curvatures 

o p 

TjW ^ V / \ 

and , the actual locations of the forces (planes I and II) must 

dx 2 dy 2 

he known if constants (such as E x , n' x , etc.) which depend upon the 
locations of the applied forces are to he evaluated. 


EVALUATION OF a AND p 
Experimental Evaluation 


The coefficients a, P, a', and p' occurring in the equations 
for the elastic constants express the effectiveness of a rih for resisting 
deformations other than bending and stretching in its longitudinal direc- 
tion. For the evaluation of cc and p for a given set of rihs 
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(longitudinal, transverse, or skew) probably sufficient accuracy will 
te achieved from a direct experimental measurement with a simple model 
having one set of ribs whose cross section and spacing duplicate those 
of the ribs for which the coefficients a and p are being sought and 
with a value of tg equal to that of the actual plate. 

A double specimen of the type shown on the right-hand side of fig- 
ure 5 may first be used to evaluate p through a tension test and, then, 
one-half of the specimen may be used to evaluate a through a bending 
test, as illustrated on the left-hand side of figure 5 . The use of a 
double specimen for the stretching test is suggested because the symmetry 
will eliminate localized bending of the skin between ribs and facilitate 
the measurement of over-all strain. Because of the prevention of local- 
ized bending, the value of p should be somewhat higher than that which 
would be obtained by stretching a single specimen like the one on the 
left-hand side of figure 5 . However, such an overestimate of p may be 
desirable if the actual plate has ribs in more than one direction, because 
then the localized curvatures associated with one set of ribs will tend 
to be reduced by the presence of the other ribs. 


The length-to-width ratio of the specimen should be great enough 
so that any end grips or heavy end sections will offer negligible resis- 
tance to transverse contraction in the stretching test and to the devel- 
opment of transverse curvature in the bending test. Furthermore the width 
of the specimen should be sufficiently large compared to the rib spacing 
so that the percentage of the specimen subject to shear-lag effects 
arising at the rib ends is small. 

The use of these tests for the evaluation of a and p will now 
be described in detail. For ease in discussion, the ribs whose a and 
p are being sought will be assumed to be oriented in the y-direction as 
shown in figure 5. After the values of Oy and Py have been deter- 
mined, however, the subscript y should be changed to x or s if , 
in the actual plate, the ribs under consideration are oriented in the 
longitudinal or skew direction of the plate. 


The conditions of the stretching test illustrated in the right-hand 
side of figure 5 are = = N, r = 0. Substituting these conditions 

a* 2 ay 2 y 

in equation (10) and making use of equations ( 31 ) , (4l), (43), (44), 
and (45) gives 
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A\j j^>x 

where, for rectangular rihs with circular fillets, — — is as given 
by equation ( 56 ) . 

By using for N x / EHe x in the right-hand side of this equation the 

value obtained in the stretching test, an experimental value of Py, 
or a, , is obtained (e x is the x-wise strain averaged over at least 
^exp 

one multiple of b x ) . 

The conditions of the bending test illustrated in the left-hand side 
of figure 5 are N x — Ny = My = 0. Substituting these conditions in 

equation (l) and making use of equations (13)> (l9)> (^ 2 ), (*+7); (^9), 
( 51 ), and ( 53 ) gives 


d£w = _ *<x 
dx 2 D x 
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( 69 ) 
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Substituting for M x /eH^ — ^ the value obtained, in the bending 


test, and for Py the value obtained from equation (68) permits equa- 

2 W 

tion (71) to yield an experimental value of cn, ( is the x-wise 

Vdx 2 

curvature averaged over at least one multiple of b x The quantities 
A g 2 , A x , Ay, A s , ky, Iy, I s are obtained from equations (4l), 
(43), (44), (45), (48), (52), and (53), respectively, with 
A Wv = a w q = J W X = x w q = °» thus > 


A s 2 = A xAy - A s 2 


** + Py 


A w v \ A A W /b x 

—21 A v = — + —II 

H y p H 


A s " 


d t S 


1 - p c 


1 % l h * r. 


K = ^ — k w 

2 Kr H w y 


I v = — + 


k W v k y A £ 


’yTs 
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12(1 - d 2 )^ 


where ku is as given by equation (59) • 


Theoretical Evaluation 

Accurate theoretical analysis of the situations depicted in figure 5 
is difficult. However, it is possible to obtain values of a and p 
that underestimate or overestimate the stiffness of the specimens. 

An underestimate is obviously obtained by assuming no part of the 
rib to be effective in resisting transverse stretching or bending in a 
direction transverse to itself. A lower-limit value of p is, therefore, 

P = Pll = 0 (73) 


When (3 Is taken as zero t 


immaterial. 
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An overestimate is obtained by analyzing the two specimens shown 
in figure 5 for their small deformations under the assumption that plane 
sections perpendicular to the skin and perpendicular or parallel to the 
direction of ribbing remain plane. The results of such an analysis of 
the two situations illustrated in figure 5 are as follows: 



for the double specimen on the right-hand side of figure 5; and 



for the single specimen on the left-hand side of figure and where I, 

g, and f are geometric properties of segments of length b x of the 
cross sections shown in figure 5* The letter I represents the moment 
of inertia of such a segment about its centroid, g is the integral, 
taken in the x-direction, of the reciprocal of the local thickness 
measured in the z-direction, and f is tg 2 times a similar integral 
of the cube of the reciprocal of the local thickness. When the ribs are 
rectangular with circular fillets, these quantities are given by the 
following formulas: 
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b x 

f = -£ - 


G w, 
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tw yf b s V 

to V h/ 


+ f ' 


(78) 


where g' and f' are functions of the ratio of fillet radius to skin 
thickness plotted in figure 6, and is as given by equation (62) 


The values of N x /EHe x and mJeH 5 — obtained from equations (jk) 
, I dx 

and (75) may be thought of as experimental results and they may therefore 
be substituted in equations (68) and (71) to obtain v al ues of f3m- and 
a UL corresponding to an overestimate of the stiffness of the specimen 


A lower overestimate of stiffness can be obtained by anal yzing, on 
the basis that plane sections remain plane, the single specimen on the 

left-hand side of figure 5 for both N x /EHe x and mJeH? and thus 

. / 8x 2 

including the localized bending that occurs during stretching. Besides 
being more conservative, the resulting values of and would 

also be more appropriate if, in the actual plate under consideration, 
there were really only one set of ribs. An upper-limit analysis con- 
ducted entirely on the specimen on the left-hand side of figure 5 would 
yield the following expression to be used in place of equation (74): 




b s/H 

1 h®\% 

4 f / b x 


+ 



(79) 


where h is tg times the integral, taken over a length b x in the 
x-direction, of the square of the reciprocal of the local thickness: 
for circular-filleted rectangular-section ribbing. 
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*x _ _ 

ts t s 


o r Wy 

2 — ■£- + 

■ts 



+ h' 


(80) 


where h' 
for Mx/eh 


is plotted in figure 

3 i£w. 

^ 2 


6 . 


Equation (75) would still he used 


EVALUATION OF a' AND p’ 


The coefficients a' and p' , which define the effectiveness of 
a rih in resisting twisting and shearing relative to its longitudinal 
and transverse directions, are not as readily measured experimentally 
nor as readily hounded hy an upper limit as a and p, although, of 
course, a lower-limit stiffness is obtained hy equating p' to zero. 

An approximate evaluation of a' and p' may he made hy assuming 
that the same volume of rih material resists shear as resists transverse 
stretching, that is, 


P'y = Py (81) 

and then hy computing where this material must he placed (a') in order 
to give the proper torsional stiffness as determined with the aid of 
reference 6. The computation of a' will now he described in detail. 


Consider an element, like the one on the left-hand side of figure 5> 
having only y-wise ribbing and subjected to a pure Mxy loading. From 
equations (3), (15) and (54) one can solve for a' v in terms of the 

I A 2 W 

measured or computed ratio j as follows: 
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(82) 
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where 


(3 


, V 


k xy ' t 
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i a 'y ) 


s . B , flw y/ t * 

T + p if~ 


(83) 


Solving for a'y gives 


a 



hi 1 , h\ 

“xy i(V \ 5 

y ~ 

/ 2 L, Aw y/ bx ts 

EH^ d 2 w 6\h / 
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X k H / 

2(1+ |a) dx by 


(84) 


The value of the ratio M* 


EH- 


8 2 


w 


to be inserted in the 


2(l + n) 8x 8y 

above formula can, in the absence of test data, be derived by an adapta- 
tion of the method used in reference 6 for computing the torsional stiff- 
ness of I-beams and H-beams, which gives 



2(l + n) bx by 



where d is the diameter of the largest circle which can be inscribed 
in the cross section at the junction of the rib and skin and can be com- 
puted from the formula 
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and a is a constant whose value depends on t w ^tg and ry^tg. The 

value of a is obtainable from figure 7 of reference 6 or, whenever 
t w /r w \ 

— * y ^ 0.6l - 0.23 — i , from the following formula: 


a = 0.09^ + 0.070 t — ~ 


The meanings of the various terms within the parentheses of equa- 
tion (85) are apparent: represents the contribution of the skin, 

considered as an infinite plate, to the twisting stiffness of the waffle; 

ifl _ ^S.\ (— 1 is similarly representative of the twisting 

5 v /tw u 


stiffness of the rib; the term with -0.105\^y corrects for the fact 

that the rib is actually not infinitely deep; and the term with 

represents the additional stiffness due to the fillets. The value 0.105 

2b w 

is based on the assumption that — -^ = 2.3; for values of 2by /ty 

t W y *1 * 

less than 2.3, the number 0.105 should be replaced by the number obtain- 
able in figure 3 of reference 6 with the abscissa label b/ n replaced 
by the label 2b^ty y . 

COMPARISON OF CALCULATED AND EXPERIMENTALLY MEASURED 
VALUES OF ELASTIC CONSTANTS 


As a partial check on the theory, experimental measurements were made 
of the stretching stiffness Eq, bending stiffness Dx, shearing stiff- 
ness Gk, and twisting stiffness Dxy of plates with integral ribs running 
either longitudinally or transversely (fig. 1 (a)) or skewed (fig. 1 (c)). 

The procedures used for the measurement of Dx and Dxy were essentially 
the same as those described in reference 5 for sandwich plates. The meas- 
urements of Eq and Gk were made with long-gage -length resistance-type 

wire strain gages mounted in the four corners, or diagonally on the four 
sides, of square-tube compression or torsion specimens similar to the 
square tubes of reference 4 . The compression specimens were tested in the 
1,200,000-pound-capacity testL^MO^^^ne^id the torsion specimens in the 
combined load testing machin4fl||^^HMB structures research laboratory. 
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The experimental values obtained, for the stiffnesses are indicated 
by the circles in figures 7 and 8. In figure 7 the stiffnesses are 
plotted against the angle of skew of the ribbing (with 0=0° and 
9 = 90 ° corresponding to purely longitudinal and purely transverse 
ribbing, respectively) for plates having nominally the same weight. In 
figure 8, for a given angle of skew (0 = 45°), the variation of the 
elastic constants with skin thickness is plotted. The relatively large 
scatter in the test data is due to the fact that the plates used were 
sand castings and, hence, had appreciable variations in thicknesses from 
one specimen to another and also within each specimen. 

For comparison, theoretical values of the four elastic constants 
were computed from equations (3l)> (l3)> (33)> and (15) and. are plotted 
in figures 7 and 8. The lowest curve in each graph is obtained from 
the lower-limit assumption, p = 0; the highest curve gives calculated 
upper-limit values based on the use of equations (74) and ( 75 ) in calcu- 
lating 07 ^ and p^; the middle (dashed) curve shows the results 

obtainable by using for a and p values determined experimentally on 
specimens like those in figure 5 * In each case it was assumed that 
p' = p, and a' was computed from equations (84) and ( 85 ). Table I 
summarizes the upper-limit and experimental values of a and p used 
for these calculations. 

In general, figures 7 and 8 indicate that the agreement between 
calculation and experiment is within the experimental scatter, with the 
calculations based on the values a exp and p exp giving the best results. 

CONCLUDING REMARKS 


On the basis of an idealization of integrally stiffened plates to 
more uniform plates resembling plywood, formulas have been derived for 
the elastic constants of- the plates with integral ribbing in one or more 
directions. Two sets of elastic -constant formulas have been given, based 
on two different forms of the force -distortion equations. 

The formulas for the elastic constants involve four coefficients a, 
P, a', and p' for each rib which define the effectiveness of the rib 
in resisting stretching and bending in its transverse direction, hori- 
zontal shearing, and twisting. Experimental means of determining these 
coefficients are discussed, as are theoretical methods of obtaining 
values corresponding to lower-limit or upper-limit assumptions regarding 
the stiffness of the plate. 
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The predictions of the formulas for four of the elastic constants 
are compared with experiment and good correlation is obtained when 
experimentally determined values (or, in most cases, upper-limit values) 
of a and p are used in the formulas for the elastic constants. 
Despite experimental scatter, the calculations and experiments agree, 
in general, both in magnitude and in regard to trends resulting from 
variation in angle of skew of ribbing or in skin thickness. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va. , May 26 , 1955 - 
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APPENDIX A 
NOMENCLATURE 


Plane I is defined as the plane in which N x acts and in which e 
is measured. Plane II is defined as the plane in which IJ y acts and iS 
Which Ey is measured. Plane III is defined as the plane in which 
acts and in which 7^ is measured. 


'XX 


xy 

“-yx 

: yy J 



c 22 J 


General Symbols 


coupling elastic constants associated with bending and 
stretching and defined by the force -distortion equa- 
tions (1), (2 ) , ( 4 ), and (5), lb -1 q 


coupling elastic constants associated with bending and 
stretching and defined by the force-distortion equa- 
tions (7), (8), (10), and (11), in. 


k coupling elastic constant associated with twist and 

shear and defined by the force-distortion equa- 
tions (9) and (12), in. 

v X ’ tending stiffnesses in x- and y-directions, respec- 

D l> d 2 j tively, in-lb * 

E xy> D k twisting stiffnesses relative to x- and y-directions 

i n — 1 "h 


Young's modulus of material, psi 


E 

E 


x> 

1 ' 



extensional stiffnesses in x- and y-directions 
tively, lb /in. 7 


re spec - 


G 


xy 


shear modulus of material, psi 

shear stiffness of plate in xy-plane, lb/in. 
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M x , M y 
N x 


s 

t 

T 


w 

U, V 
x 

y 

z 


7 


xy 


resultant bending -moment intensity in x- and 
y-directions, respectively, lb 

resultant twisting-moment intensity with, regard to 
x- and y-directions, lb 

intensity of resultant normal force acting in 
x-direction in plane I, lb/in. 

intensity of resultant normal force acting in 
x-direction in plane II, lb/in. 

intensity of resultant shear force acting in x- and 
y-directions in plane III, lb/in. 

coordinate, measured parallel to skewed rib, in. 

coordinate, measured perpendicular to skewed rib, in. 

coupling elastic constant associated with twist and 
shear and defined by the force-distortion equa- 
tions (3) and (6), lb“-*- 

displacement in z-direction, in. 

strain energy, in- lb 

coordinate, measured in longitudinal direction, in. 

coordinate, measured In transverse direction, in. 

coordinate, measured perpendicular to faces of skin, in. 

shear strain, with respect to x- and y-directions, of 
plane III 


e x> e y 


strain of plane I In x-direction and of plane II in 
y-direction, respectively 


H 


Poisson's ratio for material 


Px; Hy 


^'x> ^'yl 
Hi, M-2 J 


Poisson's ratios associated with bending in x- and 
y-directions, respectively, and defined by the force- 
distortion equations (l), \ 2 ) , ( 7 ), and (8) 

Poisson's ratios associated with extension in x- and 
y-directions, respectively, and defined by the force- 
distortion equations (4), ( 5 ), (10), and (ll) 
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^X' "by 

b s 

b S 


b w 

d 

h, k 

H 

r ¥ 

t 

0 


Symbols Representing Dimensions 

x-wise and y-wise length, respectively, of smallest 
repeating unit of plate, in. 

spacing of skew ribs, equal to b x /sin 0 or byjcos 0, 
in. 

rib spacing (measured between center lines of parallel 
ribs), in. 

rib depth, H - tg, in. 

diameter of largest circle that can be inscribed in 
cross section at intersection of rib and skin, in. 

distance from planes of zero strain to rib centroids, 
in. 

over-all height of rib plus skin, in. 

radius of fillet, in. 

thickness, in. 

angle of skewed ribbing, measured from the longitudinal 
direction, deg 


Symbols Used in Equations for Elastic Constants 


a 


constant used in equations for calculating a'jjg 


A W X ' A V Aw s 


"W 


cross-sectional area (including fillets) of x-wise, 
y-wise, and skewed ribs ^A^ includes area of two 

ribs^, sq in. 


general symbol for Ay x , 


\> or A w s 


f, g, h constants used in equations for calculating 0^ 

and Sul 

cross-sectional moment of inertia of x-wise, y-wise, 
or skewed ribs about their centroids (l w is twice 

\ s ], 

the m ome nt of inertia of a single skew rib>, in. 


I Wx > ^y* T W S 
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kjj, k m 
^ ^Wg 



dimensionless distance from middle surface of sheet 
to planes I, II, and III, respectively, expressed as 
fractions of the over-all height H 


dimensionless distance from middle surface of sheet 
to centroid of x-vlse, y-wise, or skewed rib, 
expressed as a fraction of the over-all height H 



constants used to locate the effective centroid of a 
rib for resisting bending in its transverse direction 


a 

a' x , a' y , a' s 
a' 

Pll> ^expJ ^ul 
Px> Py> Ps 


general symbol representing a x , oty, or a s 

constants used to locate the effective centroid of a 
rib for resisting twisting 

general symbol representing a' x , a 1 y , or ot 1 Q 

constants used to define effectiveness of a rib in 
resisting stretching in its transverse direction 


P general symbol representing Py, or p s 

p 1 , p* } p 1 constants used to define effectiveness of a rib in 

X ^ s resisting shearing 

P r general symbol representing p J x , p T y, or p T g 


Subscripts 

L longitudinal 

S sheet or skin 

T transverse 

W rib (web) 

s y x, y indicate application to skewed, x-wise, or y-wise ribs 

or directions 


LL 

lower limit 

UL 

upper limit 

exp 

experimental 
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E2 = 




1 - 


Ey 


*y 1 - 


c^y/ 


: xy + My c yy) + Cyy^ _ ^ y j( C yy + ^ C xy) 


Er 


E = 

y 


1 + Eo < Ci 


L 

C 12 " ^ C 22 

+ C22 

C 22 " ^ C 12 1 

1 c 12 

Eil 1 " M x^yl 

D2(i - HxVty) J 


Gfc - 


J xy 


1 - 2D X yG X yT^ 


G xy _ 


Gk 1 ^ 


°k + C k\ 


> 


(B5) 

(b6) 


D. 


C,, - C. 


x \ . „ / Ey 


'll ■ xxl L _ 


"XX 


VxVy/ 

C 11 - ^ C 21 
D l ( 1 - Wy) 


+ HxC 


, / “y 

*v - ^ 


> 


(B7) 


c i2 = C ^(i . ^ y j + ^ C yy^i - 


El 


hc^Vy 


c xy _ 


C]_2 - ^x c 22 

D l( 1 " hc^) 


(B8) 


D, 


C 21 = ^y c xx\ 


c 22 = Vy C xy\Y 


1 - 


+ C. 


Hx^yy 


y* ! - 


Wy) 


®x \ ( ^y 




°yy\r 


^y 


_ c 2i - ^y c n 

^ = M 1 - >w) 

_ c 22 ~ ^y c 12 

^ “ M 1 - L>v) 


(B9) 


(BIO) 
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APPENDIX C 

DERIVATION OF FORMULAS FOR ELASTIC CONSTANTS 


The basic assumptions of the analysis have already been described. 
In the derivations that, follow, where the word "rib" is used, it will 
usually be understood to mean one of the substitute sheets, depending 
on which property of the rib is under consideration. Separate deriva- 
tions are given for the constants associated with bending and stretching 
arid those associated with twisting and shear. 


Constants Associated With Bending and Stretching 


In the derivation of the formulas for the elastic constants asso- 
ciated with bending and stretching, an element of the integrally stiffened 
plate will be considered; the element has the average prescribed curva- 


d 2 w d 2 w / 

tures and and the strains e x (measured in some arbitrary 


c)x £ 


dy 2 


plane which will be referred to as plane I) and e y (measured in some 


other arbitrary plane which will be referred to as plane II). The 
development of these prescribed deformations requires the application 
of moments of intensity M x and My and forces of intensity N x 

(acting in plane I) and Ny (acting in plane II). These moments and 
forces and the locations of planes I and II are shown in figure 9* 


If the strains are assumed to vary linearly through the thickness 

d 2 w 


of the element, two horizontal planes can be found (in terms of 


d 2 w 


dx 2 ’ 


: X' 


dy' 


2’ 


and 


in which the x-wise strain and y-wise strain, respec- 


tively, are zero. These planes are Indicated in figure 10. 


Strains of components of plate .- The longitudinal extensional strains 
of the ribs measured at their cross-sectional centroids can be written in 
terms of the curvatures and the distance between the rib centroids and 
the planes of zero extensional strains. The strains of the x-wise, 
y-wise, and skewed ribs are, respectively. 


e 


W- 


X L 




(Cl) 
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(C2) 


% T = h l C ° s2e + k l Sln20 
S L dx^ by* 


where the subscript L denotes longitudinal direction of a rib; the 
subscript x, the x-wise rib; the subscript y, the y-wise rib; and the 
subscript s the skew rib. The distances h^, kj, h-^, and k-^ are 

shown in figure 10. 

The transverse strains of the ribs are as follows: 

£ % - -(*2 - «*> 


% T * -( h2 ' (C5) 


£ = -( h 2 - a s H )^| s1 ” 29 - (*2 

S T \ ox V 


a s H 


8^ 

8y : 


cos^9 


(C 6) 


The extensional strains of the sheet midplane 
tures are 


in terms of the curva- 


e s 


x 



(C7) 


e 



-k 2 


b^ w 

dy 2 



(C8) 
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The curvatures 

tures of the x-vise 
skew ribs is 


and 

3x 2 Sy 2 

and y-wise ribs. 


of the element are also the curva- 
respectively. The curvature of the 


^ ^ cos 2 0 + & sin 2 0 

8s 2 8x 2 8y 2 


The horizontal shear strain in one of the skew ribs, relative to 
the longitudinal and transverse directions of the rib, can be written in 
terms of the x-wise and y-wise strains at the same level, which in turn 
are determined by the x-wise and y-wise curvatures; thus. 






“' S H) 


dy 2 


sin 9 cos 9 


The x-wise and y-wise ribs have no shear strain. 


(CIO) 


Expressions for the dimensions h^, hg, h^ , k-^, k 2 , and k^.- 

In the derivation of equations (Cl) to (C8) and of equation (CIO), the 
assumption was made that the strains varied linearly from the planes of 
zero strain. On the basis of the same assumption, expressions are 
written for the strains in planes I and II - the planes in which N x 
Ny act and in which e x and e-y are measured. These expressions are 




(cn) 




(C12) 


ho = k T H - -2L_ 

* 1 a 2 w 

8x2 


(C13) 



from which 
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kg = k U H 

. fy_ 
a 2 w 

(C14) 


8y 2 


dimensions h-j_. 

hj, k-L, and 

kj may be written 

h l = k w s H 

- hg 

(C15) 

> 

\M 

ll 

- *2 

( ci6) 

k l “ k W s H 

- kg 

(C17) 

ii 

*r' 

- kg 

(Cl8) 


where k w H, k w H, k w H locate the centroidal axes of the ribs from 

the center line of the sheet. Substituting for h 2 and k 2 from equa- 
tions (C 13 ) and (CllO gives 


= ( k W s " k l) 


k I H + 


ax 2 


h 3 _ ( k w x “ k ij H + 


a g w 

ax 2 


■ ( k w s ■ 


*11) H + 


Sy 2 


k* = 


- (s • k ^) H + & 


8y c 


(C19) 


(C20) 


(C21) 


(C22) 
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Evaluation of strain energy .- The total strain energy of the element 
of the integrally stiffened plate can he written as the sum of the strain 
energies of its component parts; thus, 


V = i 


e w 2 EAtj dx + - 
w x L "x 2 



e w 2 EA w dy + i 

w y L w y 2 


fb x sec9 


0 




-EA Wg ds + 


1 

2 




+ | 



1 

2 


Pb x sec9 


0 


£ W ST 2EPsAw s dS 


+ i 


ph x sec9 

'0 


7 Ws G P's A W s dS + 


1 

2 




fs x 2 + 6 s/ + 2 “ e Sx E S y )]T7 tsdXdy + 


1 

2 



/ 5 2 w 

\dx 2 / " x 


+ 



'c> 2 v \ 2 

hj EI *s is + 



+ Ittf + 2ll 


kdx 


idy 2 / 


d 2 w d 2 w 
dx 2 dy 2 


E 


1 - p‘ 


tq3 

5F 


dxdy 


(C23) 


In equation (C23) the first three terms give the energy of exten- 
sion of the ribs in their longitudinal directions, the second three 
terms the energy of extension of the ribs in their transverse direc- 
tions, the seventh term the energy associated with the shearing of the 
ribs, and the eighth term the energy of extension of the skin. The 
next three terms give the energy of bending of the ribs, and the final 
term gives the energy of bending of the skin. 

Carrying out the integrations of equation (C23), dividing by b x by 

to reduce the result to strain energy per unit area, and substituting 
the previously derived expressions for the distortions £u > j and 

x y 

so forth gives 
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V 

^x^y 



1 — t s + + Pv + ^^(cos k 0 + p sin^e + 

. „2 S *y ^ *>x ‘bsV 5 


1- H' 

B' c — - — sin 2 0 cos 2 0] 
s 1 + u > 


e 2 + 2 
x 


— - t c + -^^sin 2 0 cos 2 0 + 


p S' v, 
1 - / Ds 


B c sin 2 0 cos 2 0 - P'o — - — sin 2 0 cos 2 ©) 
Ps s 1 + u / 


£x5y + 2 irr 7 


t s( k I H ) 


A w ./- 


— ”l k w " 

V w x I 


A Wv / \ A Wc; [7— 

H + Py -^(“y " k i) H + "71 [r w s " 


kj H cos^Q + 


H s (a s - kj)H sin^e ♦ P' s (a' s - sin 2 0 co S 2 e)]U x 0 


€ v w " + 

dx 2 


-M- 


Li - ^ 


k S ( k II^0 + ^[^Ws - k Il) H sin2e cos2q + 


Ps( a s ~ k n) H sin20 cos 2 0 - 
P's(«'s - k n) H (rTT; sln2e cos S 


a 2 


V 


1_ tg ♦ e, ^ + a. 

2 b ^ by bx b s \ 


1 - p 


cos^0 + 


p' s — - — sin 2 0 cos 2 0 

1 + H 


e y 2 + 2 { ts ' kTU 1 + 

vj- - n 


!( k I H ) 


A W S /- 


(k Wg - kj^H sin 2 0 cos 2 0 + P s (a s - k^H sin 2 0 cos 2 0 - 


tion continued on next page) 
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e's( a 's - k i) H (rr^ sin?e cos20 )} e y H + 2 {r^ ts ( kllH ) + 

Px ^( a x - k Il) H + "T^( kw y " k Il) H + ^ [*. - k n) H sin4e + 

p s (a s - k n) H cos4e + p's(a's - k n) H (rri: sin2e cos20 ]| e y + 

< f - r - — t g ^ + cos 4 0 + i_— Ukjl) 2 + 

\l 2 (l - n 2 ) S ^ *s 1 - M 2 bV 1 ' 

( k W x - k x ) 2 H 2 + Py ^(Oy - kj) 2 ^ + ^[[k Ws - kj) 2 H 2 cos^ 0 + 
p s (a s - kj^ 2 H 2 sin k 9 + P' s (a' s - k l) 2 H 2 (rf“^ sin20 cos20 ] 
y tg^ + sin 2 0 cos 2 0 + ^ 2 tg^kjkjjH 2 ) 


Aw x /r 



ll 2 (l - U 2 
^w R IT— 


~ ( k w s " k l)( k w s “ k ll) f^sin 2 © COS 2 0 + 


P s ( a s - k l)( a s - kn)H 2 sln 2 0 cos 2 0 - 

P' *l)(«'s - *nK(i7Y sm 2 e cos 2 e)]|00 

{ / 1 -~gT t s 5 + ^ + ^ sin4e + pr t s( k n H ) 2 + 

(12^1 - x D s 1 - V- 



42 



MCA RM L53E13a 


where the identities sin 0 and i = r-^- cos 9 have been sub- 

b x bg by b s 

stituted to simplify the expressions. 

Invoking the principle of virtual displacements by differentiating 
the energy expression (C24) with respect to each of the strains and 

curvatures and dividing by EH or EH 2 gives the following expressions 
for the forces and moments: 


dV' _1_ _ Nx 
de x EH EH 


_L_ *6 ♦ ^ + . t NAUs-e ♦ p s sin“e ♦ 

-L _ 2 H H ^ H H \ s 


S' — § — sin 2 0 cos 2 0J 
5 1 + u / 


e x + 


: s Aw s/ bs / 2 2 

b ■ -(sin^0 cos^ + 


2 H H 


, 2 a _2 


sin^ cos^9 - P'g — - — sin 2 9 cos 2 0 

s s 1 + p 


€y+ 1rr7 H k i + 




COS^0 + 


Ps (o^s - sin^e + P’ s (a' s - kj )[ - 2 ~ sin 2 9 cos 2 0j !>H + 


dx 2 


A W S /^ 


■±-ik TTt " B/ ^f(k w -k TT ) 
n H LA Ws J1 ) 


sin 2 0 cos 2 0 + 


P s (a s - kjj^ sin 2 9 cos 2 9 - P ’ s (a ’ s - k TI )^ sin 2 0 cos 2 0) j>H 


(C25) 
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dV _1_ = JV 

^€y EH EH 


M- tg 




(sin 2 9 cos 2 6 


+ 


P s sin 2 9 


cos 2 9 - 


P' s — - — sin 2 9 cos 2 6) 

1 + |i / 



A W x /^>y 

P* + 

^ H 




P s cos^G + p'g ■ - sin 2 G cos 2 oj 


€ y + 



AWg / 1 


H 


■E% - 4 ) 


sin 2 9 cos 2 9 + 




[•jsin^G cos 2 G - P'g^ct'g - — sln 2 9 cos^Q^H 


U S 

H 


-=£ *n + ftc 


A W x/V 
H 







k 

sin G +. 




cos^G 



s _ kj-rW — 2 — sin 2 G cos 2 g) >H ^ 

Ilj U^ ij ay 2 


(C26) 
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kk 

av i 

cPv eh 2 

3 



5c. 

EH £ 


rrH ¥ k i + ^Nc - k i) + ^ - k i) + 


A «s/ b = 17c 


H 


- kj^ cos^0 + 3 S ( a s “ k i) sin^Q 


-Tl 


k I + 


p 's( a 's - k i)(dhr sin20 cos 20 J} ex + {rr^ it 

~TL ~ |K ” k l) sin2 9 cos 2 0 + Ps^ctg - kjj sin 2 0 cos 2 0 


p ’s(“'s - k i)(r+“s sln2e oos2e |} £ y + {^d 1 . ^i Ctt) 3 + 

+ cos^e + -i-g ^ k/ + ^(E Wx - k z ) 2 

h5 h5 1 - M- n n A 


Aw y/ D x / v 2 A w s / b 6 rr_ » 

Py ii (°y - k i) + 5 [O^Ws ' k i) 

P s ( a s - kj)2 sin^0 + 3's(a' s - sin 2 0 cos 2 ©) j>H 0 + 


l2 COS^0 + 


+ Ws ^ S sin 2 0 cos 2 0 + — — ^ — k T k-rx + 


12 (l - n 2 )^ H>/ H3 1 - n 2 H 1 11 


A Ws /^ 


H 


( k W s - k l)( k W s - knjsln 2 © cos 2 0 + 


Ps(a s - k l)( a s " k n)sin 2 0 cos 2 0 - 

P's(a's - k l)( a 's - k n)(r4n sln2e cob26 |} h 0 < C2 T) 
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1 - n- 


2 Y' k n + ' Ws / ' h " k ii) sin2e cos2q + 


Ps( a s _ kjj]sin 2 0 cos 2 0 - P' s (a' s - — sin 2 0 cos 2 0jjje x + 


-1 ts 


1 - H' 


2l k II + p x — 5 — v^ x “ "II 


A Wx /by 


- k TT ) + 


A Wv/ b x,- 


r~( k w y - k u) + 


~1T ~ |fr s ~ k Il) sln4e + Ps(«s - kn)co S 4 9 

p' s (a’s - k n)(rri; s1 " 20 COB20 3} e y + 


12(1 - ^2) 


— W — — S sin 2g cos 20 + 


1 - u' 


— ICyk-rj + 

2 H I II 


W ^ S ( k W s - k l)( k W s - kjjsin 2 © cos 2 0 + 

Ps(«s - k l)( a s - k n )sin 2 0 cos 2 0 - 

p’»(“'. - k i)( a ’s - k n)(rT7 sl “ 20 c ° s20 )} H + 

/—l ft? + + fHsft sin k 9 ♦ _i_ £(„_)< 

12(l - |i 2 ) ' H / h 3 g3 1 _ ^2 H \ II/ 
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The equations for N x , Ny, M x , and My (eqs. (C25) to (C28) ) can 
he written as 


— = A x e x 
EH X X 


+ A s e y + A x (k x - k^H ^ + A s ( k s - kjjjH |-g 


^ = a s € x + Vy + A s( k s - k i) H H + VK - k n) H H (C50) 

ox ay 


“^2 “ Ax( k x " k i) € x + A s( k s " k l) e y + [^x + Axfcx-" k i)^| H 


| s + A s (k s - kj)(k B - knj H 


(C31) 


M _ 

= A s( k s " k Il) e x + Ay( k y - k Il) € y + 


+ A a (k s - kj)(k, - ^ + [ly + Ay{*y - 


where A x , A, and so forth, are given in equations (k3) to (5 k ). 


(C32) 


In order to identify the desired elastic constants associated with 
extension and bending, the foregoing force-distortion relationships, 
equations (C29) to (C32), need only to be put into the form of equa- 
tions (1), (2), (10, and (5) or (7), (8), (10), and (ll). 


Constants Associated With Twisting and Shearing 

The derivation of the formulas for the elastic constants associated 
with twisting and shearing is a p arallel one to that for the lending and 
stretching constants, 
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An element of the integrally stiffened plate will he considered 

-s2 

which has the average prescribed twist - - and shear strain 7^ 

dx dy ^ 

(measured in some arbitrary plane which will be referred to as plane III). 
These prescribed deformations can be effected by the application of 
twisting moments of intensity and shearing forces of intensity Nxy 

(acting in plane III) to the element (see fig. 11 ). 

If the horizontal shear strain is assumed to vary linearly through 

( d 2 w 

the thickness, the horizontal plane can be found I in terms of 

V dx dy 

and ?xy) h 3 ’ 8 zero shear strain. This plane is shown in figure 12 . 


Strains of components of plate .- The extensional strains of the 
longitudinal and transverse and one of the skew ribs in their longitu- 
dinal directions at their centroids are 


€ W =0 
Wx L 

(C35) 

e W = 0 
yL 

(C34) 

ey = ±h' ^w sin 20 

Ws L 1 dx dy 

(C35) 

The transverse strains of the ribs are 


£ Wy = 0 

*T 

(C36) 

O 

II 

(C37) 

e w = +( h* - a s H)-^- sin 20 

w s T V 2 s Jbx Sy 

(C38) 
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The extenslonal strains of the sheet are 


e 


Sx 


= 0 



a 2 * 

The twist g - 
ax dy 

one of these ribs 


causes bending of the diagonal ribs, 
is given by 


(C39) 

(C40) 

The curvature of 


& = - sin 20 (C4l) 

as 2 5x 

The curvatures of the longitudinal and transverse ribs are zero. The 
shear strain in the skin middle surface is given by 


7 


S 


= -2h' 2 


-.2 

a w 

ax ay 


(C42) 


The magnitude of the shear strain of the diagonal ribs is given by 


7 «s 




a 2 w 
ax ay 


cos 20 


The shear strain of the x-wise and y-wlse ribs is given by 


(C43) 


' 7 „ x * - 2 ( h '2 
7w y ’ -4 h '2 

Expressions for the dimensions 

expressions can be written for the 
which Nxy acts and in which y^ 


- 

(C44) 

- a' H) 

(C45) 

y 'ax ay 

h'^ and h^-- 

The following 


trains in plane III (the plane in 
is measured, see fig. 11): 


?xy 



- k m H ) 


a^ 

ax ay 


(C46) 



from which 


7M 


MCA RM L53E13a 



b9 


By geometry 


h'l = k W s H - h* 2 


(C48) 


Substituting for h^ from equation (C47) gives 


h ’l " ( k W s " k IIl) 


k TTX \ H + 1 


7xy 


2 

chc ^y 


(C49) 


Evaluation of strain energy .- The total strain energy can be 
written as 


u -i 


fb x sec8 


K) 


e u 2 EA U ds + i 
w s L w s 2 


'b x sec0 


0 


%,> A «s 4s + 


1 1 K T&r** + 1 j, ' + 




1 

2 


'b x sec9 rb x fby 

y W s G P’ B %s dS + | j[j J 0 ’’S G V* dy + 


0 


0 J 


'' >y G is! ^ + 

\bx by 6 


0 


1 r b * sec8 /£aV 

2 Jo \a s 2 / 


EXrj ds 
s 


(C50) 


In equation (C50) the first term gives the energy of extension of 
the skewed ribs in their longitudinal directions, the second term the 
energy of extension of the skewed ribs in their transverse directions, 
the next three terms the energy of shearing of the ribs, and the sixth 
term the energy of shearing of the skin. The next term represents the 
energy of twisting of the skin, and the last term gives the energy of 
bending of the skew ribs. 

Carrying out the integrations of equation (C50), dividing by b x by, 
substituting previously derived expressions, and so forth, gives 



ro|t*J 
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— — = U 1 
b x t) y 


+ 


1 , Aw x 

. T P'x + 


1 p-„ ^ + 


2(l + n) S " 2(1 + n) F x by ' 2(1 + n) K y b x 


— — [sin 2 0 cos 2 0 + p s sln 2 0 cos 2 0 + P'o — — r cos 2 

v 1 8 8 2(1 + ]l) 


bs L_ 


cos 20 


7xy + 


4 ( 2 d"i -) ts ( k m H ) + 2(1 + n) p ’ x ^( a ' x ■ k m) 
2(1 7 n) P V “^( a 'y " klII ) H + _ bf{( k Ws " ^n) 


H + 


H sin 2 0 cos 2 0 + 


^s( a s “ k m)® s i n ^® cos 2 0 + 


»■«(«■. - k m) H 


2(1 + \i) 


cos 2 20 


y jfci + f 1 +5 + 

7x y dy \6(1 + y) S 


Iff, 


k -bf sin2e cos 6 + rnr t s( k m H ) + 


2 o, 

1 + n P x b. 


Aw, 


( a ’x " k ni ) 2H2 + 377-37 p, y “^(“’y " ^n ) 202 + 

IT 

l ^( kw s ~ k ui) 2 H 2 sin 2 0 cos 2 0 + P S ^Q! S - kjxxj 2 H 2 sin 2 0 cos 2 0 + 


S T s (a’ 


s " k III 


N2h2 F — 1 cos 2 20 

) [2(1 + \i) 



d 2 ** 


JJ / \dx Sy, 




(C 51 ) 
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Differentiating the energy expression (C51) with respect to each 

of the distortions and dividing by EH or EH 2 gives the following 
expressions for the forces and moments: 

dU' 1 _ N xy 
d 7xy EH EH 


1 “S , 1 p. Vbr , 1 W** , 

2(1 + p) H 2(1 + p) P x H 2(1 + p) y H 


A wJ*s| 


sin 2 0 cos 2 9 + (3 S sin^G cos^G + 




cos 2 2g] 

(1 + H) J 


V + 2 2(1 + TO •# kl11 + 


1 p. V b * 

2(1 + p) P x H 


(a' x - 


l p t % /" A / a , \ 

tut p y “h r y " k m) 


“ krTTi + 


A W s / b s 


^k w - kjj^sin 2 0 cos 2 0 + 0 8 (o: s - kjjj^sin 2 © cos^G + 


coA !lf E 5T^ 
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, ^ E 112 
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2 ^SL 
EH 2 


= 2 


-1 


-§• k TTT + 


i A W x A>y , 

1 n f i fry r 


12(1 + n) T “HI - 2(1 + n) p ' x H 


( a ’x - k IIl) + 


1 Q I N/^, . X , 

iirrm p y ”h— v a y - k m) + 


A W s /*£ 

H 


■j(k Ws - k Ii;[ )sin 2 0 cos 2 0 + P B (a s - k IZI ) sin 2 0 cos 2 0 + 


2(1 + *i) 


cos^O 


P's( a 's - k IIl) 

4 ■ Ws / — sin 2 9 COS 2 0 + — ^ — ■^•( k m) 2 + 
h 3 1 + n H \ J ~ L - L ' 

2 ^x/V/ \2 + 

— P X H \ * ‘ 

2 p, A “y i> rTT \2 + 

TTZ e y y 1I1 - 1 


l7xy + 1 6(1+ ToCh) + 


%L 


4 -2£L — -!(&w s * k lll) 2sln20 cos2e + 


P s ^a s - k I3;I ) 2 sin 2 9 cos 2 0 + 


M°'s - k m) 2 [|(rT^) cos22 ! 


H 


d 2 w 
dx by 


(C53) 
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The equations for and Myy (eqs. (C52) and (C53)) can be 


written as 


N *y _ a ~ . oa Hr v \ H ^ 2 w 


" A xy?xy + 2A xy( k ^y " k m) ] 


6x dy 


(C54) 


2 eP ~ 2Ax y( kxy " k m) 7 xy + [jxy + ' k IIl)^] H ^° 55 ^ 

where A^y, k^, and I^y. are given in equations (46), ( 50 ), and (54), 
respectively. 

Equations (C54) and (C55) may readily be put into the form of 
equations (6) and (3) or (12) and (9) to yield either the original or 
the new elastic constants, respectively. 
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TABLE I 

VALUES OF a, a', p, AND p' USED IN THE CALCULATION OF 
THE ELASTIC CONSTANTS FOR COMPARISON WITH EXPERIMENTAL 
MEASUREMENTS OF E x , D x , AND D^ 
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b w /t s 

a exp 

a UL 

P f -Pexp 

“V-P| ul 

Pexp 

Pul 

'bw/t'S = °- 2 
(a) 

1 

2 

4 

8 


0.25 

.15 

.085 

.046 


0.25 

.24 

.51 

0.20 

.23 

.14 

.12 

0.53 

.45 

.29 

.19 

^wAs = 0>1 * 
00 

1 

2 
4 
8 


0.14 


0.44 


0.l4 


^These values, computed from equations ( 68 ), ( 7 l)> ( 7 ^- ) j ( 75 )> 
(84), and ( 85 ) were used for calculating constants for all configura- 
tions given in figures 7 and 8 except those for which 9=0° and 
0 = 90 ° (one-way stiffening). 

4’rhese values, computed from equations ( 68 ), (7l)> (75)> (79)> 
(84), and ( 85 ) were used for calculating constants for configurations 
of figure 7 having 9=0° and 9 = 90° • 
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(a) Longitudinal or transverse 


(b) Longitudinal and transverse 



(c) Skewed (d) Skewed plus longitudinal and transverse 


Figure 1.- Ribbing configurations considered. 


Plane I 

Plane E 

Plane M 


Figure 2.- Forces and moments acting on element. 
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Figure 4.- Comparison of idealized and actual rib-skin combinations. 





Figure 6.- Values of f ' , g' , and h' used in upper-limit evaluation 

of a and 0 . 
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having E = 10.7 x 10^ ksi, \i = 0-32, and having the following pro- 
portions: — = 4, — = 2, — = 2; in addition, for 0 * 0° or 90°, 
t s t w ry 

^ = 0.4 and for 0° < 9 < 90°, SS ■ 0.2. 
t»S b S 




= 0<2> ^ = 0.2 in., E = 10.7 X 10 3 ksi, n = 0.52. 

b S 
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Plane I 

Plane IT 

Figure 9 *- Forces and moments considered for analysis of bending and 

stretching. 


L 




Plane of zero x-wlse strain 

Plane of zero y-wise strain 

® Rib centroids 

Figure 10.- Dimensions for analysis of "bending and stretching. 
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Plane ffl 


Figure 11.- Shears and moments considered for analysis of twisting and 

shearing. 
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Plane of zero shear strain 

® Rib centroids 

Figure 12.- Dimensions for analysis of twisting and shearing. 
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